Mapping the functional impact of non-coding regulatory elements in primary T cells through single-cell CRISPR screens.

BACKGROUND: Drug targets with genetic evidence are expected to increase clinical success by at least twofold. Yet, translating disease-associated genetic variants into functional knowledge remains a fundamental challenge of drug discovery. A key issue is that the vast majority of complex disease associations cannot be cleanly mapped to a gene. Immune disease-associated variants are enriched within regulatory elements found in T-cell-specific open chromatin regions. RESULTS: To identify genes and molecular programs modulated by these regulatory elements, we develop a CRISPRi-based single-cell functional screening approach in primary human T cells. Our pipeline enables the interrogation of transcriptomic changes induced by the perturbation of regulatory elements at scale. We first optimize an efficient CRISPRi protocol in primary CD4+ T cells via CROPseq vectors. Subsequently, we perform a screen targeting 45 non-coding regulatory elements and 35 transcription start sites and profile approximately 250,000 T -cell single-cell transcriptomes. We develop a bespoke analytical pipeline for element-to-gene (E2G) mapping and demonstrate that our method can identify both previously annotated and novel E2G links. Lastly, we integrate genetic association data for immune-related traits and demonstrate how our platform can aid in the identification of effector genes for GWAS loci. CONCLUSIONS: We describe "primary T cell crisprQTL" - a scalable, single-cell functional genomics approach for mapping regulatory elements to genes in primary human T cells. We show how this framework can facilitate the interrogation of immune disease GWAS hits and propose that the combination of experimental and QTL-based techniques is likely to address the variant-to-function problem.

Use of Bulk Segregant Analysis for Determining the Genetic Basis of Azole Resistance in the Opportunistic Pathogen Aspergillus fumigatus.

A sexual cycle was described in 2009 for the opportunistic fungal pathogen Aspergillus fumigatus, opening up for the first time the possibility of using techniques reliant on sexual crossing for genetic analysis. The present study was undertaken to evaluate whether the technique 'bulk segregant analysis' (BSA), which involves detection of differences between pools of progeny varying in a particular trait, could be applied in conjunction with next-generation sequencing to investigate the underlying basis of monogenic traits in A. fumigatus. Resistance to the azole antifungal itraconazole was chosen as a model, with a dedicated bioinformatic pipeline developed to allow identification of SNPs that differed between the resistant progeny pool and resistant parent compared to the sensitive progeny pool and parent. A clinical isolate exhibiting monogenic resistance to itraconazole of unknown basis was crossed to a sensitive parent and F1 progeny used in BSA. In addition, the use of backcrossing and increasing the number in progeny pools was evaluated as ways to enhance the efficiency of BSA. Use of F1 pools of 40 progeny led to the identification of 123 candidate genes with SNPs distributed over several contigs when aligned to an A1163 reference genome. Successive rounds of backcrossing enhanced the ability to identify specific genes and a genomic region, with BSA of progeny (using 40 per pool) from a third backcross identifying 46 genes with SNPs, and BSA of progeny from a sixth backcross identifying 20 genes with SNPs in a single 292 kb region of the genome. The use of an increased number of 80 progeny per pool also increased the resolution of BSA, with 29 genes demonstrating SNPs between the different sensitive and resistant groupings detected using progeny from just the second backcross with the majority of variants located on the same 292 kb region. Further bioinformatic analysis of the 292 kb region identified the presence of a cyp51A gene variant resulting in a methionine to lysine (M220K) change in the CYP51A protein, which was concluded to be the causal basis of the observed resistance to itraconazole. The future use of BSA in genetic analysis of A. fumigatus is discussed.

Modifying the m6A brain methylome by ALKBH5-mediated demethylation: a new contender for synaptic tagging.

Synaptic plasticity processes, which underlie learning and memory formation, require RNA to be translated local to synapses. The synaptic tagging hypothesis has previously been proposed to explain how mRNAs are available at specific activated synapses. However how RNA is regulated, and which transcripts are silenced or processed as part of the tagging process is still unknown. Modification of RNA by N6-methyladenosine (m6A/m) influences the cellular fate of mRNA. Here, by advanced microscopy, we showed that m6A demethylation by the eraser protein ALKBH5 occurs at active synaptic ribosomes and at synapses during short term plasticity. We demonstrated that at activated glutamatergic post-synaptic sites, both the YTHDF1 and YTHDF3 reader and the ALKBH5 eraser proteins increase in co-localisation to m6A-modified RNAs; but only the readers showed high co-localisation to modified RNAs during late-stage plasticity. The YTHDF1 and YTHFDF3 readers also exhibited differential roles during synaptic maturation suggesting that temporal and subcellular abundance may determine specific function. m6A-sequencing of human parahippocampus brain tissue revealed distinct white and grey matter m6A methylome profiles indicating that cellular context is a fundamental factor dictating regulated pathways. However, in both neuronal and glial cell-rich tissue, m6A effector proteins are themselves modified and m6A epitranscriptional and posttranslational modification processes coregulate protein cascades. We hypothesise that the availability m6A effector protein machinery in conjunction with RNA modification, may be important in the formation of condensed synaptic nanodomain assemblies through liquid-liquid phase separation. Our findings support that m6A demethylation by ALKBH5 is an intrinsic component of the synaptic tagging hypothesis and a molecular switch which leads to alterations in the RNA methylome, synaptic dysfunction and potentially reversible disease states.

TCRs with Distinct Specificity Profiles Use Different Binding Modes to Engage an Identical Peptide-HLA Complex.

The molecular rules driving TCR cross-reactivity are poorly understood and, consequently, it is unclear the extent to which TCRs targeting the same Ag recognize the same off-target peptides. We determined TCR-peptide-HLA crystal structures and, using a single-chain peptide-HLA phage library, we generated peptide specificity profiles for three newly identified human TCRs specific for the cancer testis Ag NY-ESO-1157-165-HLA-A2. Two TCRs engaged the same central peptide feature, although were more permissive at peripheral peptide positions and, accordingly, possessed partially overlapping peptide specificity profiles. The third TCR engaged a flipped peptide conformation, leading to the recognition of off-target peptides sharing little similarity with the cognate peptide. These data show that TCRs specific for a cognate peptide recognize discrete peptide repertoires and reconciles how an individual's limited TCR repertoire following negative selection in the thymus is able to recognize a vastly larger antigenic pool.

Gfi1aa and Gfi1b set the pace for primitive erythroblast differentiation from hemangioblasts in the zebrafish embryo.

The transcriptional repressors Gfi1(a) and Gfi1b are epigenetic regulators with unique and overlapping roles in hematopoiesis. In different contexts, Gfi1 and Gfi1b restrict or promote cell proliferation, prevent apoptosis, influence cell fate decisions, and are essential for terminal differentiation. Here, we show in primitive red blood cells (prRBCs) that they can also set the pace for cellular differentiation. In zebrafish, prRBCs express 2 of 3 zebrafish Gfi1/1b paralogs, Gfi1aa and Gfi1b. The recently identified zebrafish gfi1aa gene trap allele qmc551 drives erythroid green fluorescent protein (GFP) instead of Gfi1aa expression, yet homozygous carriers have normal prRBCs. prRBCs display a maturation defect only after splice morpholino-mediated knockdown of Gfi1b in gfi1aa qmc551 homozygous embryos. To study the transcriptome of the Gfi1aa/1b double-depleted cells, we performed an RNA-Seq experiment on GFP-positive prRBCs sorted from 20-hour-old embryos that were heterozygous or homozygous for gfi1aa qmc551 , as well as wt or morphant for gfi1b We subsequently confirmed and extended these data in whole-mount in situ hybridization experiments on newly generated single- and double-mutant embryos. Combined, the data showed that in the absence of Gfi1aa, the synchronously developing prRBCs were delayed in activating late erythroid differentiation, as they struggled to suppress early erythroid and endothelial transcription programs. The latter highlighted the bipotent nature of the progenitors from which prRBCs arise. In the absence of Gfi1aa, Gfi1b promoted erythroid differentiation as stepwise loss of wt gfi1b copies progressively delayed Gfi1aa-depleted prRBCs even further, showing that Gfi1aa and Gfi1b together set the pace for prRBC differentiation from hemangioblasts.

Nanopore sequencing and assembly of a human genome with ultra-long reads.

We report the sequencing and assembly of a reference genome for the human GM12878 Utah/Ceph cell line using the MinION (Oxford Nanopore Technologies) nanopore sequencer. 91.2 Gb of sequence data, representing ∼30× theoretical coverage, were produced. Reference-based alignment enabled detection of large structural variants and epigenetic modifications. De novo assembly of nanopore reads alone yielded a contiguous assembly (NG50 ∼3 Mb). We developed a protocol to generate ultra-long reads (N50 > 100 kb, read lengths up to 882 kb). Incorporating an additional 5× coverage of these ultra-long reads more than doubled the assembly contiguity (NG50 ∼6.4 Mb). The final assembled genome was 2,867 million bases in size, covering 85.8% of the reference. Assembly accuracy, after incorporating complementary short-read sequencing data, exceeded 99.8%. Ultra-long reads enabled assembly and phasing of the 4-Mb major histocompatibility complex (MHC) locus in its entirety, measurement of telomere repeat length, and closure of gaps in the reference human genome assembly GRCh38.

MinION Analysis and Reference Consortium: Phase 2 data release and analysis of R9.0 chemistry.

BACKGROUND: Long-read sequencing is rapidly evolving and reshaping the suite of opportunities for genomic analysis. For the MinION in particular, as both the platform and chemistry develop, the user community requires reference data to set performance expectations and maximally exploit third-generation sequencing. We performed an analysis of MinION data derived from whole genome sequencing of Escherichiacoli K-12 using the R9.0 chemistry, comparing the results with the older R7.3 chemistry. METHODS: We computed the error-rate estimates for insertions, deletions, and mismatches in MinION reads. RESULTS: Run-time characteristics of the flow cell and run scripts for R9.0 were similar to those observed for R7.3 chemistry, but with an 8-fold increase in bases per second (from 30 bps in R7.3 and SQK-MAP005 library preparation, to 250 bps in R9.0) processed by individual nanopores, and less drop-off in yield over time. The 2-dimensional ("2D") N50 read length was unchanged from the prior chemistry. Using the proportion of alignable reads as a measure of base-call accuracy, 99.9% of "pass" template reads from 1-dimensional ("1D")  experiments were mappable and ~97% from 2D experiments. The median identity of reads was ~89% for 1D and ~94% for 2D experiments. The total error rate (miscall + insertion + deletion ) decreased for 2D "pass" reads from 9.1% in R7.3 to 7.5% in R9.0 and for template "pass" reads from 26.7% in R7.3 to 14.5% in R9.0. CONCLUSIONS: These Phase 2 MinION experiments serve as a baseline by providing estimates for read quality, throughput, and mappability. The datasets further enable the development of bioinformatic tools tailored to the new R9.0 chemistry and the design of novel biological applications for this technology. ABBREVIATIONS: K: thousand, Kb: kilobase (one thousand base pairs), M: million, Mb: megabase (one million base pairs), Gb: gigabase (one billion base pairs).

Signatures of Selection for Environmental Adaptation and Zebu × Taurine Hybrid Fitness in East African Shorthorn Zebu.

The East African Shorthorn Zebu (EASZ) cattle are ancient hybrid between Asian zebu × African taurine cattle preferred by local farmers due to their adaptability to the African environment. The genetic controls of these adaptabilities are not clearly understood yet. Here, we genotyped 92 EASZ samples from Kenya (KEASZ) with more than 770,000 SNPs and sequenced the genome of a pool of 10 KEASZ. We observe an even admixed autosomal zebu × taurine genomic structure in the population. A total of 101 and 165 candidate regions of positive selection, based on genome-wide SNP analyses (meta-SS, Rsb, iHS, and ΔAF) and pooled heterozygosity (Hp) full genome sequence analysis, are identified, in which 35 regions are shared between them. A total of 142 functional variants, one novel, have been detected within these regions, in which 30 and 26 were classified as of zebu and African taurine origins, respectively. High density genome-wide SNP analysis of zebu × taurine admixed cattle populations from Uganda and Nigeria show that 25 of these regions are shared between KEASZ and Uganda cattle, and seven regions are shared across the KEASZ, Uganda, and Nigeria cattle. The identification of common candidate regions allows us to fine map 18 regions. These regions intersect with genes and QTL associated with reproduction and environmental stress (e.g., immunity and heat stress) suggesting that the genome of the zebu × taurine admixed cattle has been uniquely selected to maximize hybrid fitness both in terms of reproduction and survivability.

Development of Gene-Based SSR Markers in Winged Bean (Psophocarpus tetragonolobus (L.) DC.) for Diversity Assessment.

Winged bean (Psophocarpus tetragonolobus) is an herbaceous multipurpose legume grown in hot and humid countries as a pulse, vegetable (leaves and pods), or root tuber crop depending on local consumption preferences. In addition to its different nutrient-rich edible parts which could contribute to food and nutritional security, it is an efficient nitrogen fixer as a component of sustainable agricultural systems. Generating genetic resources and improved lines would help to accelerate the breeding improvement of this crop, as the lack of improved cultivars adapted to specific environments has been one of the limitations preventing wider use. A transcriptomic de novo assembly was constructed from four tissues: leaf, root, pod, and reproductive tissues from Malaysian accessions, comprising of 198,554 contigs with a N50 of 1462 bp. Of these, 138,958 (70.0%) could be annotated. Among 9682 genic simple sequence repeat (SSR) motifs identified (excluding monomer repeats), trinucleotide-repeats were the most abundant (4855), followed by di-nucleotide (4500) repeats. A total of 18 SSR markers targeting di- and tri-nucleotide repeats have been validated as polymorphic markers based on an initial assessment of nine genotypes originated from five countries. A cluster analysis revealed provisional clusters among this limited, yet diverse selection of germplasm. The developed assembly and validated genic SSRs in this study provide a foundation for a better understanding of the plant breeding system for the genetic improvement of winged bean.

Transcriptomic responses of mixed cultures of ascomycete fungi to lignocellulose using dual RNA-seq reveal inter-species antagonism and limited beneficial effects on CAZyme expression.

Gaining new knowledge through fungal monoculture responses to lignocellulose is a widely used approach that can lead to better cocktails for lignocellulose saccharification (the enzymatic release of sugars which are subsequently used to make biofuels). However, responses in lignocellulose mixed cultures are rarely studied in the same detail even though in nature fungi often degrade lignocellulose as mixed communities. Using a dual RNA-seq approach, we describe the first study of the transcriptional responses of wild-type strains of Aspergillus niger, Trichoderma reesei and Penicillium chrysogenum in two and three mixed species shake-flask cultures with wheat straw. Based on quantification of species-specific rRNA, a set of conditions was identified where mixed cultures could be sampled so as to obtain sufficient RNA-seq reads for analysis from each species. The number of differentially-expressed genes varied from a couple of thousand to fewer than one hundred. The proportion of carbohydrate active enzyme (CAZy) encoding transcripts was lower in the majority of the mixed cultures compared to the respective straw monocultures. A small subset of P. chrysogenum CAZy genes showed five to ten-fold significantly increased transcript abundance in a two-species mixed culture with T. reesei. However, a substantial number of T. reesei CAZy transcripts showed reduced abundance in mixed cultures. The highly induced genes in mixed cultures indicated that fungal antagonism was a major part of the mixed cultures. In line with this, secondary metabolite producing gene clusters showed increased transcript abundance in mixed cultures and also mixed cultures with T. reesei led to a decrease in the mycelial biomass of A. niger. Significantly higher monomeric sugar release from straw was only measured using a minority of the mixed culture filtrates and there was no overall improvement. This study demonstrates fungal interaction with changes in transcripts, enzyme activities and biomass in the mixed cultures and whilst there were minor beneficial effects for CAZy transcripts and activities, the competitive interaction between T. reesei and the other fungi was the most prominent feature of this study.

Real-time selective sequencing using nanopore technology.

The Oxford Nanopore Technologies MinION sequencer enables the selection of specific DNA molecules for sequencing by reversing the driving voltage across individual nanopores. To directly select molecules for sequencing, we used dynamic time warping to match reads to reference sequences. We demonstrate our open-source Read Until software in real-time selective sequencing of regions within small genomes, individual amplicon enrichment and normalization of an amplicon set.

Characterisation of the horse transcriptome from immunologically active tissues.

The immune system of the horse has not been well studied, despite the fact that the horse displays several features such as sensitivity to bacterial lipopolysaccharide that make them in many ways a more suitable model of some human disorders than the current rodent models. The difficulty of working with large animal models has however limited characterisation of gene expression in the horse immune system with current annotations for the equine genome restricted to predictions from other mammals and the few described horse proteins. This paper outlines sequencing of 184 million transcriptome short reads from immunologically active tissues of three horses including the genome reference "Twilight". In a comparison with the Ensembl horse genome annotation, we found 8,763 potentially novel isoforms.

The role of carbon starvation in the induction of enzymes that degrade plant-derived carbohydrates in Aspergillus niger.

Fungi are an important source of enzymes for saccharification of plant polysaccharides and production of biofuels. Understanding of the regulation and induction of expression of genes encoding these enzymes is still incomplete. To explore the induction mechanism, we analysed the response of the industrially important fungus Aspergillus niger to wheat straw, with a focus on events occurring shortly after exposure to the substrate. RNA sequencing showed that the transcriptional response after 6h of exposure to wheat straw was very different from the response at 24h of exposure to the same substrate. For example, less than half of the genes encoding carbohydrate active enzymes that were induced after 24h of exposure to wheat straw, were also induced after 6h exposure. Importantly, over a third of the genes induced after 6h of exposure to wheat straw were also induced during 6h of carbon starvation, indicating that carbon starvation is probably an important factor in the early response to wheat straw. The up-regulation of the expression of a high number of genes encoding CAZymes that are active on plant-derived carbohydrates during early carbon starvation suggests that these enzymes could be involved in a scouting role during starvation, releasing inducing sugars from complex plant polysaccharides. We show, using proteomics, that carbon-starved cultures indeed release CAZymes with predicted activity on plant polysaccharides. Analysis of the enzymatic activity and the reaction products, indicates that these proteins are enzymes that can degrade various plant polysaccharides to generate both known, as well as potentially new, inducers of CAZymes.

The dynamics of genome replication using deep sequencing.

Eukaryotic genomes are replicated from multiple DNA replication origins. We present complementary deep sequencing approaches to measure origin location and activity in Saccharomyces cerevisiae. Measuring the increase in DNA copy number during a synchronous S-phase allowed the precise determination of genome replication. To map origin locations, replication forks were stalled close to their initiation sites; therefore, copy number enrichment was limited to origins. Replication timing profiles were generated from asynchronous cultures using fluorescence-activated cell sorting. Applying this technique we show that the replication profiles of haploid and diploid cells are indistinguishable, indicating that both cell types use the same cohort of origins with the same activities. Finally, increasing sequencing depth allowed the direct measure of replication dynamics from an exponentially growing culture. This is the first time this approach, called marker frequency analysis, has been successfully applied to a eukaryote. These data provide a high-resolution resource and methodological framework for studying genome biology.

Accelerated growth in the absence of DNA replication origins.

DNA replication initiates at defined sites called origins, which serve as binding sites for initiator proteins that recruit the replicative machinery. Origins differ in number and structure across the three domains of life and their properties determine the dynamics of chromosome replication. Bacteria and some archaea replicate from single origins, whereas most archaea and all eukaryotes replicate using multiple origins. Initiation mechanisms that rely on homologous recombination operate in some viruses. Here we show that such mechanisms also operate in archaea. We use deep sequencing to study replication in Haloferax volcanii and identify four chromosomal origins of differing activity. Deletion of individual origins results in perturbed replication dynamics and reduced growth. However, a strain lacking all origins has no apparent defects and grows significantly faster than wild type. Origin-less cells initiate replication at dispersed sites rather than at discrete origins and have an absolute requirement for the recombinase RadA, unlike strains lacking individual origins. Our results demonstrate that homologous recombination alone can efficiently initiate the replication of an entire cellular genome. This raises the question of what purpose replication origins serve and why they have evolved.

Genome-wide transcriptional response of Trichoderma reesei to lignocellulose using RNA sequencing and comparison with Aspergillus niger.

BACKGROUND: A major part of second generation biofuel production is the enzymatic saccharification of lignocellulosic biomass into fermentable sugars. Many fungi produce enzymes that can saccarify lignocellulose and cocktails from several fungi, including well-studied species such as Trichoderma reesei and Aspergillus niger, are available commercially for this process. Such commercially-available enzyme cocktails are not necessarily representative of the array of enzymes used by the fungi themselves when faced with a complex lignocellulosic material. The global induction of genes in response to exposure of T. reesei to wheat straw was explored using RNA-seq and compared to published RNA-seq data and model of how A. niger senses and responds to wheat straw. RESULTS: In T. reesei, levels of transcript that encode known and predicted cell-wall degrading enzymes were very high after 24h exposure to straw (approximately 13% of the total mRNA) but were less than recorded in A. niger (approximately 19% of the total mRNA). Closer analysis revealed that enzymes from the same glycoside hydrolase families but different carbohydrate esterase and polysaccharide lyase families were up-regulated in both organisms. Accessory proteins which have been hypothesised to possibly have a role in enhancing carbohydrate deconstruction in A. niger were also uncovered in T. reesei and categories of enzymes induced were in general similar to those in A. niger. Similarly to A. niger, antisense transcripts are present in T. reesei and their expression is regulated by the growth condition. CONCLUSIONS: T. reesei uses a similar array of enzymes, for the deconstruction of a solid lignocellulosic substrate, to A. niger. This suggests a conserved strategy towards lignocellulose degradation in both saprobic fungi. This study provides a basis for further analysis and characterisation of genes shown to be highly induced in the presence of a lignocellulosic substrate. The data will help to elucidate the mechanism of solid substrate recognition and subsequent degradation by T. reesei and provide information which could prove useful for efficient production of second generation biofuels.

Planarian MBD2/3 is required for adult stem cell pluripotency independently of DNA methylation.

Planarian adult stem cells (pASCs) or neoblasts represent an ideal system to study the evolution of stem cells and pluripotency as they underpin an unrivaled capacity for regeneration. We wish to understand the control of differentiation and pluripotency in pASCs and to understand how conserved, convergent or divergent these mechanisms are across the Bilateria. Here we show the planarian methyl-CpG Binding Domain 2/3 (mbd2/3) gene is required for pASC differentiation during regeneration and tissue homeostasis. The genome does not have detectable levels of 5-methylcytosine (5(m)C) and we find no role for a potential DNA methylase. We conclude that MBD proteins may have had an ancient role in broadly controlling animal stem cell pluripotency, but that DNA methylation is not involved in planarian stem cell differentiation.

Transient exposure to low levels of insecticide affects metabolic networks of honeybee larvae.

The survival of a species depends on its capacity to adjust to changing environmental conditions, and new stressors. Such new, anthropogenic stressors include the neonicotinoid class of crop-protecting agents, which have been implicated in the population declines of pollinating insects, including honeybees (Apis mellifera). The low-dose effects of these compounds on larval development and physiological responses have remained largely unknown. Over a period of 15 days, we provided syrup tainted with low levels (2 µg/L(-1)) of the neonicotinoid insecticide imidacloprid to beehives located in the field. We measured transcript levels by RNA sequencing and established lipid profiles using liquid chromatography coupled with mass spectrometry from worker-bee larvae of imidacloprid-exposed (IE) and unexposed, control (C) hives. Within a catalogue of 300 differentially expressed transcripts in larvae from IE hives, we detect significant enrichment of genes functioning in lipid-carbohydrate-mitochondrial metabolic networks. Myc-involved transcriptional response to exposure of this neonicotinoid is indicated by overrepresentation of E-box elements in the promoter regions of genes with altered expression. RNA levels for a cluster of genes encoding detoxifying P450 enzymes are elevated, with coordinated downregulation of genes in glycolytic and sugar-metabolising pathways. Expression of the environmentally responsive Hsp90 gene is also reduced, suggesting diminished buffering and stability of the developmental program. The multifaceted, physiological response described here may be of importance to our general understanding of pollinator health. Muscles, for instance, work at high glycolytic rates and flight performance could be impacted should low levels of this evolutionarily novel stressor likewise induce downregulation of energy metabolising genes in adult pollinators.

High quality de novo sequencing and assembly of the Saccharomyces arboricolus genome.

BACKGROUND: Comparative genomics is a formidable tool to identify functional elements throughout a genome. In the past ten years, studies in the budding yeast Saccharomyces cerevisiae and a set of closely related species have been instrumental in showing the benefit of analyzing patterns of sequence conservation. Increasing the number of closely related genome sequences makes the comparative genomics approach more powerful and accurate. RESULTS: Here, we report the genome sequence and analysis of Saccharomyces arboricolus, a yeast species recently isolated in China, that is closely related to S. cerevisiae. We obtained high quality de novo sequence and assemblies using a combination of next generation sequencing technologies, established the phylogenetic position of this species and considered its phenotypic profile under multiple environmental conditions in the light of its gene content and phylogeny. CONCLUSIONS: We suggest that the genome of S. arboricolus will be useful in future comparative genomics analysis of the Saccharomyces sensu stricto yeasts.

Characterisation of retroviruses in the horse genome and their transcriptional activity via transcriptome sequencing.

The recently released draft horse genome is incompletely characterised in terms of its repetitive element profile. This paper presents characterisation of the endogenous retrovirus (ERVs) of the horse genome based on a data-mining strategy using murine leukaemia virus proteins as queries. 978 ERV gene sequences were identified. Sequences were identified from the gamma, epsilon and betaretrovirus genera. At least one full length gammaretroviral locus was identified, though the gammaretroviral sequences are very degenerate. Using these data the RNA expression of these ERVs were derived from RNA transcriptome data from a variety of equine tissues. Unlike the well studied human and murine ERVs there do not appear to be particular phylogenetic groups of equine ERVs that are more transcriptionally active. Using this novel approach provided a more technically feasible method to characterise ERV expression than previous studies.